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Abstract: We investigate dark matter (DM) in the context of the minimal supersym- 
metric extension of the standard model (MSSM). We scan through the MSSM parameter 
space and search for solutions that (a) are consistent with the Higgs discovery and other 
collider searches; (b) satisfy the flavor constraints from B physics; (c) give a DM candidate 
with the correct thermal relic density; and (d) are allowed by the DM direct detection 
experiments. For the surviving models with our parameter scan, we find the following fea- 
tures: (1) The DM candidate is largely a Bino-like neutralino with non-zero but less than 
20% Wino and Higgsino fractions; (2) Constraints from the Higgs sector and rare b decay 
measurements exclude the low mass region favored by the DAMA, CoGeNT, and CRESST 
experiments; (3) The relic density requirement clearly pins down the solutions from the 
Z and Higgs resonances (Z, h, H, A funnels) and co-annihilations; (4) Future direct search 
experiments will likely fully cover the Z, h funnel regions, and H , A funnel regions as well 
except for the "blind spots" ; (5) Future indirect search experiments will be more sensitive 
to the CP-odd Higgs exchange due to its s-wave nature; (6) The branching fraction for the 
SM-like Higgs decay to DM can be as high as 8%, while those from heavier Higgs decays to 
neutralinos and charginos can be as high as 20%. We show that collider searches provide 
valuable information complementary to what may be obtained from direct detections and 
astroparticle observations, and that the Higgs bosons may play an essential role. 
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1. Introduction 

Observations of the cosmic microwave background, gravitational lensing, clustering of 
galaxies, galactic rotation curves, etc. have provided compelling evidence for the existence 
of Dark Matter (DM), which is likely to be of particle origin. One of the best motivated 
candidates for DM is the Weakly Interacting Massive Particle (WIMP), a good example 
of which is the Lightest Supersymmetric Particle (LSP) (for reviews, see [1-3]). If WIMPs 
exist in the Galaxy, they may be detected through direct search experiments [4-11]. The 
DAMA experiment [4] has detected an annual modulation in the measured recoil spectrum 
at the 8.9a level, consistent with the presence of WIMP DM in the Galaxy. More recently, 
the CoGeNT [5] and CRESST [6] experiments have also obtained results that are consistent 
with low mass WIMP DM. On the other hand, these results have been challenged by other 
experiments such as XENON-10 [7], XENON-100 [8] and more recently TEXONO [12], 
which have excluded the parameter space favored by the DAMA, CoGeNT, and CRESST 
experiments. Complementary to the direct searches, indirect detection experiments in- 
clude the Fermi gamma ray space telescope [13], Air Cherenkov Telescopes [14-16], and 
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the Wilkinson Microwave Anisotropy Probe (WMAP) [17,18]. The WMAP observations 
place a lower limit on the particle mass m x > 10 GeV, assuming a velocity-independent an- 
nihilation cross section (a a v) = 1 pbxc [19-24]. The non-observation of gamma rays from 
DM annihilation in the nearby dwarf galaxies [13, 25] has been used to place constraints 
on the DM particle mass m x > 40 GeV, for neutralino annihilation to the bb channel with 
a velocity-independent cross section (cr&v) = 1 pbxc, although these bounds would be 
relaxed with a more general analysis including the velocity-dependent contributions [26]. 

On the other hand, the LHC experiments have made a historic discovery of the long- 
sought-after Higgs boson predicted by the Standard Model (SM). The experiments also 
show no evidence for Beyond-SM Higgs bosons, nor other new physics such as Supersymme- 
try (SUSY) etc. with the current data, seemingly in favor of heavy colored sparticles [27-31]. 
Several authors have recently studied the present LHC data and the implications for DM, 
as well as the possibility that future LHC data will provide information to the DM puz- 
zle [32-51]. Although the SUSY parameter space has been significantly reduced due to the 
absence of a SUSY signal at the LHC and due to the constraining properties of the SM-like 
Higgs boson, a dark matter candidate can still be readily accommodated in SUSY theories. 

With the ever increasing experimental sensitivity of DM detection experiments, we are 
motivated to explore to what extent DM properties have been constrained by the results 
from particle accelerator experiments. Our goal is to systematically examine the com- 
plementarity between DM direct detection experiments, indirect detection searches, and 
collider experiments, and in particular explore the potential pivotal role played by the 
Higgs bosons. We perform a comprehensive study in the framework of the minimal super- 
symmetric extension of the standard model (MSSM). We impose the following constraints 
on our model considerations: 

(1) Relic abundance: The neutralino LSP constitutes all the cold DM, consistent with 
the cosmological observations [17,18]. 

(2) Collider constraints: The MSSM parameter space satisfies all collider constraints 
from the Higgs boson searches and has a SM-like Higgs boson near 126 GeV. 

(3) Flavor constraints: The parameter space satisfies the flavor constraints from b —> 
S7 [52], and B s -> fi+ yT [53]. 

We further check the consistency of the annihilation rate at zero velocity {c^v) (v — > 0) with 
CMB observations, and the absence of gamma rays from nearby dwarf galaxies [13,25]. It 
is known that the spin-independent WIMP-nucleon elastic scattering cross section obtained 
by the XENON- 100 experiment [8] puts a very strong bound on the MSSM parameter space. 
We find that the surviving region has characteristic features, notably a Bino-like LSP. What 
is most interesting to us is that all these scenarios would lead to definitive predictions for 
the LHC experiments, that can be verified by the next generation of direct /indirect search 
experiments such as LUX [10] and XENON-IT [11]. 

The rest of the paper is organized as follows. In Section 2, we compute the WIMP 
relic density in a model-independent manner. We emphasize the importance of including 
the effect of velocity-dependent annihilation (a&v) = a + bv 2 + 0(t> 4 ), which is crucial 
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when the b term from the p-wave is not negligible. This will have direct consequences in 
the interpretation of indirect search results. In Section 3, we discuss our technique for 
scanning the MSSM parameter space. In Section 4, we present our results, and discuss the 
experimental constraints from the Higgs and flavor searches. We also discuss the constraints 
on the parameter space imposed by the XENON- 100 search for spin- independent scattering, 
as well as the Super-K and IceCube/DeepCore limits on spin-dependent scattering. We 
show that future experiments such as LUX and XENON-IT will likely probe the natural 
supersymmetric parametric space consistent with the LSP constituting all the DM. We 
present extensive discussions of our results in Section 5 and finally draw our conclusions in 
Section 6. Details of the relic density calculation are presented in the Appendix. 

2. Dark matter relic density 

With in the Standard Cosmology, we evaluate the thermal history of the dark matter [1,2]. 
We assume that the WIMP, generically denoted by %, constitutes all of the thermal DM. 
Let us define the DM relic density £l x as the ratio of the DM mass density at the present 
epoch to the critical mass density (/9 C rit)- The well-measured value from WMAP9, ACT, 
SPT, SNLS3, BAO, and HO [18] is obtained 

n h 2 = ra*"*, = 0.1148 ± 0.0019. (2.1) 

X Pcrit/h 2 

The number density of WIMPs at time t, n x (t), can be obtained by solving the Boltzmann 
equation (see Appendix for details). The annihilation cross section (er a ) characterizes the 
WIMP dynamics for a given theory. Since WIMPs are non-relativistic at freeze-out, the 
velocity averaged cross section (a a v) may be expanded in v, customarily written as 

- a + bv 2 + 0(v 4 ), / * 1 = a+ — with x = — f, (2.2) 



lpbxc lpbxc x T 

where the traditional units are 1 pbxc = 3 x 10~ 26 cm 3 /s. Simple threshold arguments 
indicate that s-wave annihilation contributes dominantly to a, while p-wave annihilation 
contributes only to b. Requiring that the thermal relic density satisfy the measured value 
gives us the cross section. This leads to the result: 

n x h 2 « 0.11 (cr a u) ps 2.18 x 10^ 26 cm 3 /s. (2.3) 

Fig. 1(a) shows the evolution of the number of relativistic degrees of freedom, from [54]. 
We show the WIMP number density approaching the present day value in Fig. 1(b) for 
m x = 100 GeV for the extreme cases b = and a = 0. The range of a and b values 
versus the WIMP mass is shown in Fig. 1(c). For a WIMP mass above around 4 — 5 
GeV, a remains almost constant for arbitrary WIMP masses when b is of negligible value, 
reflecting the "WIMP miracle" that leads to the correct relic density. Below 4 — 5 GeV, 
somewhat larger values of a and b would be needed to yield the correct relic density, due to 
the reduction of g at the quark-hadron transition. We will not explore the very low mass 
region any further in this work. The interplay between a and b follows a linear relation 
empirically, and is shown in Fig. 1(d) for various WIMP mass values. 
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(a) (b) 




Figure 1: (a) The number of relativistic degrees of freedom as a function of temperature. The 
relic density fl x h 2 = 0.11; (b) WIMP number density evolution with temperature for illustrative 
values of a, b with m x = 100 GeV, (c) a, b values versus the WIMP mass, (d) a versus b for m x =5, 
100, and 1000 GeV. 



3. The MSSM parameters relevant to DM studies 



In SUSY theories with conserved R-parity, the lightest supersymmetric particle (LSP) is a 
viable WIMP DM candidate. For both theoretical and observational considerations [1,2, 
55-57] , it is believed that the best candidate is the lightest Majorana mass eigenstate which 
is an admixture of the Bino (B), Wino (W3), and Higgsinos (H^ >u ), with the corresponding 
soft SUSY breaking mass parameters M%, M2, and the Higgs mixing fj,, respectively. The 
neutralino mass matrix in the Bino-Wino-Higgsino basis is given by 
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where m z is the Z boson mass, 8 W the Weinberg angle, and tan/3 = v u /vd is the ratio of 
the vacuum expectation values for the two Higgs doublets. The lightest neutralino is a 
linear combination of the superpartners 

Xi = NuB + N12W3 + Ni^Hd + NuH u , (3.1) 

where iVy are the elements of the matrix N that diagonalize M neut : 

iV*M neu tiV _1 = diag{m x o,m x o,m x o,m x o}. (3.2) 

The eigenvalues of M ncut are the masses of the four neutralinos. An interesting limit 
is m z <C I Mi ± fi\ and |M2 ± in which case, the mass eigenstates (neutralinos x^) are 
nearly pure gauge eigenstates (gauginos and Higgsinos). This also implies that large mixing 
of gaugino and Higgsino components for the mass eigenstates only takes place when M\ 
and/or M 2 are nearly degenerate with /x. We will focus only on the lightest neutralino 
(henceforth denoted by Xi) with a mass m x o. In particular, we assume that it constitutes 
the majority of the DM. 

Intimately related to the neutralinos is the Higgs sector. The tree level Higgs masses 
in the MSSM can be expressed in terms of tan/3 and the CP-odd mass Ma- Radiative 
corrections enhance the Higgs mass significantly via the top quark Yukawa coupling, the 
third generation squark mass parameters Mq 3 , M;j 3 , and the left-right squark mixing ^It- 
Flavor physics observations from the 6-quark sector often serve as stringent constraints 
and we therefore include the sbottom sector parameters Md 3 and the squark mixing A\y. 
The last potentially relevant sector is the stau, which could be light and contribute to the 
t-channel exchange, co-annihilations to control the relic density. We therefore generously 
vary the MSSM parameters in the ranges 

5GeV < |Mi| < 2000 GeV, 100 GeV < |M 2 , n\ < 2000 GeV, 

3<tan/3< 55, 80 GeV < M A < 1000 GeV, (3.3) 

-4000 GeV < A t < 4000 GeV, 100 GeV < M Q3 , M U3 < 3000 GeV, 

-4000 GeV <A h < 4000 GeV, 100 GeV < M D3 < 3000 GeV, 

-4000 GeV < A T < 4000 GeV, 100 GeV < M L3 , M E3 < 3000 GeV. 

The lowest values of Mi, M 2 and fi control the LSP mass for the WIMP DM. The lower 
values of 100 GeV for M 2 , /(/ are dictated by the LEP-2 bound from the largely model- 
independent chargino searches. The lower limit of tan j3 is close to the LEP-2 Higgs search 
exclusion. The lower limit of Ma is chosen to cover the non-decoupling Higgs sector as 
well as above the LEP-2 bound on the charged Higgs. The upper limit of Mi, M 2 , fi and 
the soft SUSY breaking masses in the stop and stau sectors are set with consideration of 
naturalness [36,58-62] . The other soft supersymmetry breaking parameters are less relevant 
for our DM considerations and we therefore set the other trilinear mass parameters to be 
zero, and the other soft SUSY breaking masses at 3 TeV. 

While the natural value of \i is supposed to be close to the electroweak scale, we vary 
fj, up to 2 TeV to capture some interesting futures such as the scenario of "well-tempered 
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neutralino" [63] . Letting \i ~ 2 TeV would already allow for a severe fine tuning at the level 
of about 0.04% [58]. Although not our focus, we have included the arbitrary signs for the 
Mi, M2, fi parameters. This allows us to see the possible solutions with very specifically 
chosen parameter relations such as the "blind spots" scenarios [48,64]. 

We choose a flat prior for the scanning with a total number of scanned points around 
10 million. Several different layers of scanning are performed to account for different 
experimental constraints, as seen by the corresponding color codes in our plots. 

4. Current Constraints and the Scanning Results 

The hints of DM detection from the DAMA, CoGeNT, and CRESST experiments have 
drawn significant interest in considering valid theoretical interpretations. The sensitivity 
of the DM direct searches have been steadily improving at an impressive pace, notably 
with the XENON collaboration [7,8]. The indirect searches from WMAP, Fermi-LAT, and 
IceCube have also played crucial roles in exploring the nature of the DM particle. 

Although the null results of searching for Supersymmetry at colliders have significantly 
tightened the viable SUSY parameter region, the bounds on WIMP DM properties are 
only limited within specific models, most notably in mSUGRA or CMSSM [32,33]. The 
direct exploration of the electroweak gaugino sector at the LHC would be very challenging 
given the hostile background environment and the current search results depend on several 
assumptions [65, 66]. On the other hand, if we demand the correct WIMP LSP relic 
abundance from the current observations as in Eq. (2.1), the SUSY parameter space of 
Eq. (3.3) will be notably constrained in the Higgs and gaugino sectors. We assume a 10% 
theoretical uncertainty in the calculation of the DM relic density [35,67]. Applying the 
combined WMAP9 + ACT + SPT + SNLS3 + BAO + HO result in Eq. (2.1) combined 
with 10% theoretical uncertainty, we demand that the relic density in our model points be 
within the following 2a window 

0.0915 < Q x oh 2 < 0.1381 (4.1) 

We use the publicly available feynhiggs code [68] as the spectrum calculator. The Higgs 
constraints are imposed using the higgsbounds package [69] with our additional modifi- 
cations. We modify the codes to include the most recent LHC constraints. 1 The standard 
SLHA [70] output recorded is then supplied to the micromegas code [71] which computes 
the DM relic density, direct/indirect search cross sections and flavor calculations. This 
is done to avoid any possible inconsistency due to the subtle differences in the spectrum 
calculator, particularly the lack of accuracy in the default approximate diagonalization 
routine for the neutralino mass matrix. 

4.1 Constraints from the Higgs searches and the flavor sector 

The discovery of a SM-like Higgs boson h as well as the upper limits on difference channels 
for the other Higgs bosons A, H, shed much light on the electroweak sector, and can 

x We have implemented the latest LHC search results from ATLAS (http://atlas.ch) and CMS 
(http://cms.web.cern.ch). 
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thus guide us for DM studies. When scanning over the SUSY parameter space as in 
Eq. (3.3), and requiring the correct WIMP LSP relic abundance to be within the 2cj 
window in Eq. (4.1), we further require the theory to have a SM-like Higgs boson, and to 
accommodate all the current constraints from the Higgs searches: 

123 GeV <m h < 128 GeV, cr 77 > 0.8 a 77 (5M), 

plus Higgs search bounds from LEP, Tevatron, LHC, (4-2) 
plus LEP bounds on the chargino mass (> 100 GeV) 
and the slepton mass (> 80 GeV). 

Due to the correlation of the Higgs couplings, the requirement of the <7 77 cross section 
effectively sets the SM-like values for ay/w-, ®zz as well. 

The absence of tree-level flavor changing neutral currents (FCNC) in the SM puts 
strong constraints on new physics. We consider two processes that have been observed to 
be consistent with the SM prediction and thus provide constraints on the MSSM parameter 
space. The first process is b — > 37 [72], for which the branching fraction is sensitive to the 
charged Higgs boson and supersymmetric particles (e.g. chargino/stop) in the loop. The 
world average of the branching fraction of this channel [73] is (3.55 ± 0.24 ± 0.09) x 10 -4 
in good agreement with the standard model prediction [74-76] (3.15 ± 0.23) x 10 -4 . 

The second process is B s — > which receives a large contribution in the MSSM 

proportional to (tan 6 [77]. The LHCb collaboration has recently announced the 

first discovery [53] of this very rare decay and the branching ratio for this process was 
found to be (3.2^}'2 t.$ 3) x 10 -9 in good agreement with the standard model prediction 
of (3.23 ± 0.27) x 10" 9 [78]. They also announced a 95% confidence level upper bound of 
5.1 x 10~ 9 on BR(i? s — > fi + fi~). We adopt this more stringent value to put constraints on 
BR(B S -> n + pL-). 

We adopt the theoretical uncertainties from the SM predictions. We note that the 
uncertainties from experiments are of the same order of magnitude as the theoretical un- 
certainty for BR(6 — > 57), and thus the latter become very important. In light of these 
precision results, we require our MSSM solutions to be within 2a of the observed value: 

2.86 x 10" 4 < BR(6 -> sj) < 4.24 x 10~ 4 , 

BR(B S /z+aO < 5.1 x 10~ 9 . (4.3) 



4.2 Confronting the direct and indirect searches 

Thus far, the most stringent constraints on the spin-independent elastic scattering cross 
section (cp 31 ) come from the XENON-100 experiment. The results from the XENON ex- 
periment challenge the signal hints from DAMA, CoGeNT and CRESST in the low mass 
region of m x ~ 10 GeV, and cut deeply into the parameter space with a^ 1 ~ 2 x 10 -9 pb 
at m x ~ 60 GeV. Limits on the spin-dependent cross section are not as constraining. We 
account for the bounds from the Super-Kamiokande [79], and the IceCubc/DcepCore [80] 
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Figure 2: Allowed parameter regions versus the CP-odd Higgs boson mass m^, for (a) tan/3, 
(b) the Higgs mixing parameter /i, (c) stop mixing parameter A t and (d) LSP DM mass m x , 
respectively. All points pass the collider and Higgs constraints of Eq. (4.2). The grey squares 
require that the DM does not overdose the Universe; the red stars in addition satisfy the flavor 
constraints of Eq. (4.3); the blue disks are consistent with the LSP being all of the DM (i.e. predicts 
the correct relic density of Eq. (4.1)). The green squares pass the XENON-100 direct search bound 
in addition to the other requirements. 

experiments that are sensitive to the spin-dependent scattering of DM with Hydrogen in 
the sun. We also take into account bounds obtained by the Fermi satellite from the absence 
of gamma rays from the nearby dwarf galaxies. 

4.3 Scanning results 

We now present our results for the allowed parameter regions in Figs. 2—9. In Fig. 2, 
we show the parameter points passing the Higgs constraints in Eq. (4.2) versus the CP- 
odd Higgs boson mass itia, for (a) tan/3, (b) Higgs mixing parameter fj,, (c) stop mixing 
parameter A t and (d) DM mass m x , respectively. These allowed parameter regions are 
shown in Fig. 3 for branching fractions (a) BR(6 — > s-y) and (b) BR(S S — > fi + fi~) versus 
tan (5. We show from the same set of points, the Wino mass parameter M<i and the Higgsino 
mass parameter ji versus the Bino mass parameter M\ in Figs. 4(a) and (b). We show the 
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Figure 3: Allowed branching fraction regions versus tan/3, for (a) b — > s-y, (b) B s — » The 
corresponding experimental central values and 2a bands are plotted on each panel. Symbols and 
legends are the same as in Fig. 2. 

second and third neutralino masses m v o, ttj v o, the light stau mass and the light stop mass 

A2 A3 

versus the LSP mass m x in Figs. 4(c), (d), (e) and (f). In the above Figures 2—4, all 
points satisfy the collider, and Higgs search requirements in Eq. 4.2. The grey squares 
show MSSM models that do not overdose the universe. The red stars in addition satisfy 
the flavor requirements in Eq. (4.3). The blue disks represent the models that give the 
correct relic density in Eq. (4.1). Finally, the green squares pass the severe XENON-100 
direct search bound on the WIMP-proton spin-independent elastic scattering. 

The results obtained here are consistent with the existing literature on the studies at 
the LHC [81,82]. We make the following important observations: 

(1). Higgs constraints (grey squares): We start with points that do not overdose the uni- 
verse and satisfy the collider search requirements in Eq. (4.2). We reproduced the known 
results that there are two surviving regions: 

(1) The non-decoupling regime where ttia ~ 95 — 130 GeV, the heavy CP-even Higgs {H ) 
is SM-like, and the light CP-even Higgs (h) is nearly degenerate in mass with the CP-odd 
Higgs (A). This region is particularly interesting since it leads to rich collider phenomenol- 
ogy and favors a light WIMP mass m x < 50 GeV. These points are not shown on the plots 
since they are disfavored by the flavor constraints, as discussed next. 

(ii) The decoupling regime where ttia > 250 GeV, the light CP-even Higgs is SM-like, and 
the heavy CP-even Higgs is nearly degenerate in mass with the CP-odd Higgs. This regime 
is difficult to observe at the LHC when uia ^ 400 GeV and tan (3 ~ 10 — 20 in traditional 
SM Higgs search channels due to severely suppressed couplings to the gauge bosons. 

(2) . Flavor constraints (red stars): The two decay processes 6—7-57 and B s — > fj, + fi~ are 
the most constraining ones. The experimental central values are plotted on the calculated 
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Figure 4: Allowed parameter regions for (a) fi and (b) M-i versus the Bino mass parameter M\. 
(c) and (d) show the second and the third neutralino masses m x o, m x o for different values of the 
DM mass, (e) and (f) show the stau mass ra? and stop mass mi respectively, versus the DM mass 
m x . Symbols and legends are the same as in Fig. 2. 



branching fractions in Fig. 3 on each panel, along with 2a bands, which is summarized 
in Eq. (4.3). These flavor constraints prefer lower tan/3 values and essentially remove the 
light Higgs (H°, A , H^) solutions in the non-decoupling region in our generic scan. The 
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solutions with a light LSP of m x < 30 GeV are also eliminated. Our results are in good 
agreement with the existing studies [38,47]. Further improvements in the B s — > /i + fj,~ 
measurement would strongly constrain the large tan/3 and low rriA region. However, we 
have not tried to exhaust parameter choices with possible cancelations among different 
SUSY contributions, and some sophisticated scanning may still find solutions with certain 
degrees of fine-tuning [83] . 

(3) . Relic density requirement (blue disks): Merely requiring that the LSP does not over- 
close the universe does not constrain the MSSM parameter space very much, as most 
clearly seen from the gray squares and red stars in Fig. 4. This is because the Higgsino-like 
or Wino-like LSPs and NLSPs can annihilate efficiently through gauge bosons and Higgs 
bosons. Requiring the correct dark matter relic density at the present epoch does constrain 
the parameter space significantly. We see the preference for fi > M\ and M2 > Mi, as in 
Figs. 4(a) and (b). Otherwise the Higgsino or Wino LSP would annihilate too efficiently, 
and result in too little dark matter. Nevertheless, we do find a nearly degenerate region 
of a Bino LSP and Wino NLSPs as seen in Fig. 4(c), which is best characterized by the 
"well-tempered" scenario [63]. This scenario, however, seems to be less implementable with 
Higgsino NLSPs as seen in Figs. 4(b) and (d), if /1 is not much greater 1 TeV. Importantly 
for our interests, we see prominent strips near m x ~ mz/2, rrih/2 which are the Z and 
Higgs funnel regions. Interestingly, as seen in Fig. 2(d), there is a region of depletion near 
rriA ~ 2m x , indicating the very (too) efficient annihilation near the ^4° funnel in the s- 
channel that is removed by the correct relic density requirement. This is a result of a lower 
bound on the LSP-Higgsino component ./V13 that we will discuss later. 

(4) . Direct search bounds (green, yellow and magenta squares): The results from DM di- 
rect searches can be translated to spin-independent cross sections and thus to the MSSM 
parameters. This is shown in Fig. 5, where all the points in the colored shaded region give 
the correct relic abundance in Eq. (4.1), satisfy the collider constraints in Eq. (4.2) and the 
flavor constraints in Eq. (4.3). The parameter space favored by the DAMA, CoGeNT and 
CRESST experiments, as well as the stringent bound from the XENON-100 experiment 
are plotted. We see that there are no surviving points in the low mass region and the blue 
region is further excluded by the XENON-100 experiment 2 . As seen in Fig. 2(a), lower 
tan/3 and higher m& values are preferred. Figures 2(b) and 4(b) show the lower bound 
H > 200 GeV. This consequently leads to a heavier X3 as seen in 4(d), while X2 could be 
still as light as the LSP x? as seen in 4(c). 

The most important observation from our study is that the surviving points are quite 
characteristic. We can identify the following classes of predictive features for the LSP DM 
from Fig. 5. 

I-A (green) X1X1 —> Z —> SM predicts m x ~ mz/2 ~ 45 GeV, the Z-funnel. 

I-B (green) X1X1 ~~ * h — > SM predicts m x ~ ~ 63 GeV, the h- funnel. 

2 It should be noted that the theoretical calculation of the spin-independent cross section may have 
significant uncertainties [84,85]. 
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Figure 5: Spin-independent cross section versus the DM mass m x o. All the points in the colored 
shaded region give the correct relic abundance in Eq. (4.1), satisfy the collider constraints in Eq. (4.2) 
and the flavor constraints in Eq. (4.3). The green region represents the model points with the Z 
and Higgs resonances. The Z funnel and h funnel regions are clearly visible for WIMP masses 
around half the Z mass and half the Higgs mass. The yellow points represent the region of co- 
annihilation with Wino-likc/Higgsino-likc NLSPs. The magenta points represent the region with 
f , £> T , b, t contributions. The gray points represent the scenarios with special cancellations when 
Mi and p take opposite signs. The DAMA and CoGeNT contours (3tr) are shown for astrophysical 
parameters vq = 220 km/s, v csc = 600 km/s, and for a local density po = 0.3 GeV/cm 3 . CRESST 
contours are 2a regions, from [6]. The blue region is excluded by the XENON-100 experiment (90% 
exclusion curve from [8], for vq = 220 km/s, v csc = 544 km/s, po = 0.3 GeV/cm 3 ). Recent results 
from the TEXONO [12] collaboration are shown. Expected exclusion bounds from the ongoing 
LUX experiment [10] and the future XENON-IT experiment [11] are also shown. 



I-C (green) XiXi -> H,A-> SM predicts m x « m A ^ H /2 ~ 0.2-0.5 TeV, the il/A-funnel. 
II- A (yellow) Neutralino/chargino coannihilation [86,87]: XiXji XiXf ~^ SM. 
II-B (magenta) Sfermion assistance [88-90]: Xi^, X?^ x\b —> SM; t-channel f, v in X^X% 

We categorize model points as scenario I if the difference between the mediator mass 
and twice the LSP mass is within 8% of the mediator mass, namely 

\m z ,h,A ~ 2m x o | < 0.08 m z , h ,A- (4.4) 
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Figure 6: (a) The annihilation cross section (a a v) in the limit v — > along with the 95% exclusion 
obtained by the Fermi satellite from the absence of gamma rays from the nearby dwarf galaxies [13]. 
(b) The spin-dependent scattering cross section with a proton, along with the 90% exclusion curves 
from the Super-K [79] and IceCube [80] experiments. Legends are the same as in Fig. 5. 

We categorize model points as scenario II-A if the difference between the LSP mass and 
neutralino NLSP 3 mass is less than 15% of the LSP mass, namely m^o — m^o < 0.15m^o. 
Other cases are categorized as scenario II-B. Our classification and categorization have 
been verified by investigating a fraction of our generated model points and looking into 
their individual contributing annihilation channels. 

(5). Indirect search bounds: 

There exist cosmological bounds from the indirect search for DM signals. We present 
the annihilation cross section (cr a t)) in the limit v — > (i.e. the f-independent component) 
versus the LSP DM mass in Fig. 6(a), along with the 95% exclusion obtained by the 
Fermi-LAT satellite from the absence of gamma rays from the nearby dwarf galaxies [13]. 
We see that further improvement from the measurement at the Fermi-LAT will reach 
the current sensitivity range. The spin-dependent scattering cross section with a proton 
is shown in Fig. 6(b), along with the 90% exclusion curves from the Super-Kamiokande 
experiment [79] and the IceCube constraint from DM annihilation in the Sun [80]. We see 
that IceCube results are cutting into the relevant parameter region closing the gap from 
the direct searches, although the bounds from the indirect searches are not quite as strong 
as that from XENON-100. 

5. Discussions 

5.1 The nature of the DM 

Experimental results from the collider searches, the 6-quark rare decays and the direct DM 
searches, combined with the relic density requirement have put very stringent constraints 

3 This is almost always true because we have a very Bino-like LSP. For cases with fi, t\ NLSP with the 
sfermion coannihilation mechanism, they fall into scenario II-B automatically. 
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Figure 7: The gaugino and Higgsino fractions of the LSP versus m x o (a) consistent with the relic 
density, collider, and flavor constraints, (b) consistent with XENON- 100 in addition to the other 
requirements. The gray points represent the results for Mi and /x to have opposite signs and the 
corresponding fractions Nfi, iV^ 2 , N% 3 and A^ 4 are denoted by hollow squares, circles, daggers and 
hollow triangles, respectively. 

on the SUSY parameter space. This in turn could have significant implications for searches 
at future collider experiments. Of primary importance is the nature of the LSP. We show 
the gaugino and Higgsino fraction (Nf { ) of the neutralino LSP versus its mass in Fig. 7, 
consistent with all collider and flavor measurements as well as the correct relic density. 
From Fig. 7(a), we note that the surviving points are mostly Bino-like (N^, as shown 
by the black dots), with lower fractions for Wino-like (Nf 2 , red dots) and Higgsino-like 
(iV^ 3 , iVj 4 , green and blue dots, respectively). As noted earlier, this is because Wino-like 
and Higgsino-like LSP's annihilate very efficiently via SU(2) gauge interactions resulting 
in too little dark matter at the present epoch. Yet, the LSP could not be purely Bino-like 
since it would overdose the Universe. In the region m v o ~ 40 GeV— 60 GeV, the line 
structures corresponding to the Higgsino components are easily identifiable with the Z and 
h exchanges. 

The XENON- 100 direct search plays a crucial role in fixing the DM properties. The 
surviving points are shown in Fig. 7(b). We see that the Wino and Higgsino fractions of the 
LSP are further constrained to a smaller fraction. The component comes in with the 
opposite signs with respect to the H u and W components. Bino-like LSP becomes more 
pronounced and the Wino and Higgsino components consist of less than 7% each, rendering 
the "well-tempered" scenario [63] difficult to realize with large Bino-Wino or Bino-Higgsino 
mixing. The comparison between Fig. 7(a) and (b) clearly shows the XENON- 100 exclu- 
sion probes deeply into the Higgsino and Wino components. On the other hand, the Hd 
component iV^ must be non-zero, and so is for H u . The non-zero lower bound would 
have significant implications for direct searches as we will discuss next, although the precise 
values may depend on the choice of the ranges for Mi, M2 and (i. 

It is important to note that a relative opposite sign between Mi and fj, could result 
in a subtle cancellation for the hxx coupling [48,64], and thus enlarge the allowed mixing 
parameters, reaching about 20% Wino/Higgsino mixtures, as shown by the grey points in 
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Fig. 7(b). This can happen only for a higher LSP mass when co-annihilations or H, A 
funnels could be in effect. 



5.2 Lower limit on the spin-independent cross section 

With our assumptions in the MSSM framework and the well-constrained properties of the 
LSP, we would expect that the DM scattering cross section may be predicted. 

It is interesting to ask whether one may derive a lower limit for the spin-independent 
scattering cross section. This is quite achievable for the Higgs resonance situation. Much 
effort has been made to derive the neutralino recoil cross sections with nuclei in various 
SUSY models [56,91-96]. This cross section mainly receives contributions from h, H 
exchanges and squark exchanges. Given the current experimental bounds on the masses of 
the squarks from the LHC [97,98], the Higgs exchanges dominate. As a good approximation 
in the decoupling limit cos(a — /3) ~ 0, tan f3 > 3 and a Bino-like LSP, the neutralino- 
nucleon cross sections via the Higgs exchanges are of the form [95] : 



It s , It u an d /tg are parameters measured from nuclear physics experiments. The H 
exchange benefits from an enhancement factor (N13/ cos/3) 2 . When the H is heavy, the h 
exchange proportional to iV 2 4 becomes important. Due to our natural choices of parameters 
as in Eq. (3.3), there exist lower bounds on and iV 2 4 , as shown in Fig. 7, although 
unnaturally large values of \x and uia could relax these bounds. Consequently, the LSP spin- 
independent cross sections at the Z, h funnels, which are mainly from the LSP scattering 
of a t-channel H exchange, reaches a lower bound, roughly 



As seen in Fig. 5, this range (green dots) will be probed by the ongoing LUX experiment 
and will be fully covered by the proposed XENON-IT experiment. Similar argument could 
be also applicable to the H, A funnel regions, where t-channel h exchange could become 
dominant. However, an exception is that a subtle cancellation takes place when M\ and 
\i take opposite signs [48,64]. As seen from the grey points in Fig. 5(b), this can happen 
only for a higher LSP mass when co-annihilations or H, A funnels could be in effect. 

In Ref. [99], a parameter-independent lower bounds a^ 1 > 2 x 10 -12 pb could be 
obtained in the mass range 440 GeV < m x < 1020 GeV and ji > 0. In the most general 
pMSSM [46] with much larger M2, n parameters, the spin-independent cross section could 
go lower, depending on the mixing parameters. 

5.3 Connection to the indirect searches 

The WIMP DM at the present epoch is non-relativistic and we can thus relate the current 
indirect search via LSP annihilation to that at freeze out [100]. The partial wave properties 
of the LSP annihilation allow us to understand the various contributions. Fig. 8(a) shows 



a X N oc < 




Ut s + §ffTG) 2 , H exchange, 
Ut u + 4ffTG) 2 , h exchange. 



(5.1) 



|JVll| 2 |JVl4| 2 



a^(h,H)>10- w pb. 



(5.2) 
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Figure 8: (a) the annihilation cross section (a a v)(v — > 0) versus the spin independent cross section 
(jp 1 . (b) the spin-dependent cross section afP versus the spin independent cross section af/ . 
Legends are the same as in Fig. 5. 
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Figure 9: Branching fractions to ncutralinos and charginos (a) for h, H, A decays to LSP pair 
versus the LSP mass, (b) for H, (c) for A, and (d) for versus its mass respectively. 



the annihilation cross section (a & v) (v — > 0) versus the spin-independent cross section when 
scattering off a proton Up . The model points in green squares near half the Z boson 



and near half the 126 GeV Higgs boson in Fig. 6(a) correspond to the low branch of 
the green squares in Fig. 8(a), due to the p-wave suppression. On the other hand, the 
s-channel annihilation through A in the mass window 200 GeV ~ 500 GeV in Fig. 6(a) is 
through s-wave, and thus has a relatively high cross section (indicated by the high branch 
of green squares). Although the LSP couplings to H and A both are mainly through 
their Higgsino components N14, the H exchange is via p-wave and thus yields a lower 
cross section as shown by the middle branch in Fig. 8(a). Finally, we note that the LSP- 
NLSP co-annihilation (yellow squares) could yield higher cross sections for both direct, and 
indirect searches, depending on their Wino and Higgsino components. Fig. 8(b) shows the 
spin-dependent cross section versus the spin-independent cross section, for our different 
models. Some of the models represented by blue points have a large enough cross section 
to be probed by IceCube/DeepCore [80]. However, a large spin-dependent cross section 
implies a proportionally large spin-independent cross section. Thus, all models in blue are 
excluded by the XENON-100 experiment. Fig. 8 illustrates the connection between spin- 
dependent and spin-independent measurements, as well as the connection between direct 
searches and indirect searches. Further improvement of the indirect search sensitivity will 
soon reach the relevant parameter region, and will probe the A exchange contribution due 
to its s-wave dominance. 



5.4 Implication of LSP for Higgs physics 

As shown in Fig. 10, a class of solutions exist with the LSP mass nearly half the mediator 
Z, h, H, A mass that undergoes a resonant enhancement in annihilation, in the Higgs 
funnel region. One may expect to see the mediator's invisible decay mode to LSP pairs 
in collider experiments. Unfortunately, these channels are kinematically suppressed near 
threshold by the non-relativistic velocity factor. Near the Z peak for example, the search 
for Z — > X1X1 would prove impossible since the branching fraction would be smaller than 
10~ 5 due to this suppression. On the other hand, invisible decay channels could be sizable 
for heavier parent particles. Shown in Fig. 9(a) are the branching fractions of h, H, A to a 
pair of LSP X1X1 versus its mass, which would be the invisible mode in collider experiments. 
It is informative to note that the SM-like Higgs boson receives two distinctive contributions 
denoted by the red squares 

RR (h _> v °v°l ~ / °- 8% m * * 60 GeV ' d SI 

BR max (h X1X1) ~ I g% ^ ^ 45 QeY (5-3) 

The branching fraction near 60 GeV is rather small although this is clearly identifiable as 
the /i-funnel region. The branching fraction near 45 GeV is about an order of magnitude 
larger because of the available kinematics, even though it is from the Z- funnel. This leads 
to the very interesting and challenging possibility of observing the Higgs invisible decay 
at the LHC [101-103], (a sensitivity of about 20% is considered feasible). The search 
sensitivity would be significantly improved at future e + e~ colliders, reaching about a few 
percent at the International Linear Collider (ILC), and even 0.3% at the TLEP [104]. 
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Table 1: Connection between the SUSY DM properties and Higgs bosons. 



5.5 Consequences of co-annihilation 

For the co-annihilation scenarios, the other lower- lying SUSY particles are nearly degener- 
ate with the LSP to ensure efficient annihilation. The common case is that the NLSP and 
NNLSP of the Winos {X\, X2) or the Higgsinos (x% > X23) are nearly degenerate with the 
Bino-like LSP, with appreciable mixing among them. On the other hand, the XENON- 100 
search bound puts a constraint on the sizes of the mixing as seen from Fig 7(a) and (b). 
Nevertheless, the spin-independent cross sections are typically higher than those from the 
Higgs resonances, reaching ~ 10 ^ pb (yellow region). The indirect detection cross 
sections are in general between the s-wave dominance (higher green band) and p-wave 
dominance (lower green band). 

As shown in Figs. 9(b)— (d), branching fractions for the other heavy Higgs bosons 
H, A to a pair of light SUSY particles could reach up to about 10% — 20%. These are 
the solutions for the correct relic density with co-annihilations. However, due to the mass 
degeneracy, the final decay products would be rather soft and are difficult to observe with 
the LHC. Consequently, these also yield the invisible decay channels. 

The coannihilation scenarios predict a rich spectrum near the LSP mass, leading to 
many different phenomena that can be explored by sparticles pair productions [66,105]. 

6. Summary and Conclusions 

Within the framework of the minimal supersymmetric extension of the standard model 
(MSSM), we investigated the possibility of the lightest supersymmetric particle being all 
the dark matter in light of the recent discovery of a SM-like Higgs boson, and the search 
for other Higgs bosons and SUSY particles at the LHC. We scanned through a wide range 
of the MSSM parameter space, and searched for model points wherein LSP has the correct 
properties to be the (WIMP) thermal DM. We studied the freeze-out of WIMPs and 
computed the cross section required to give the correct DM relic abundance in terms of 
velocity-independent and velocity-dependent components. We applied the constraints on 
the MSSM Higgs sector from the LEP, Tevatron and LHC observations. We also imposed 
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flavor constraints from the recent experimental results at the LHCb and BELLE, and found 
stringent bounds on the parameter space. The low LSP mass region is essentially closed, 
yielding a rough bound m x > 30 GeV. We thus find that the DAMA/CoGeNT/CRESST 
observations are incompatible with the MSSM. 

The XENON- 100 experiment significantly constrains the viable parameter region via 
the spin-independent elastic WIMP-proton scattering cross section. We are able to identify 
the Higgs contributions and thus to make predictions for future searches at the LHC. There 
are also clear contributions from the co-annihilation channels. Table 1 summarizes these 
distinctive MSSM model points, and the relation with the Higgs bosons. We reiterate the 
key points of our findings. For the resonance scenarios as in I-A, I-B, and I-C, 

• Z, h, H and A are the most important mediators at resonance to yield the correct 
relic abundance and give predictive narrow mass windows as shown in Fig. 5, which 
we refer to as the Z, h and if/A-funnel regions. It is interesting to note that the 
spin-independent scattering in the Z, h funnel can be dominated by the i-channel H 
exchange when A^ 3 3> N^. 

• With our parameter scanning, the necessarily non-zero Wino, Higgsino components 
of the LSP imply a lower bound for the WIMP scattering cross section mediated 
by h and H, as in Eq. (5.2). In particular, the spin-independent cross sections may 
be fully covered by the next generation of direct search experiments for DM mass 
around 30 - 800 GeV such as LUX and XENON-IT. An exception is the fine-tuned 
cancellation, the "blind spots" scenario, way above the Z, h funnels. 

• Z, h, H and A mediators determine the partial wave decomposition as listed in 
Table 1 and predict a definite range of indirect search cross sections. It is especially 
sensitive to the ^-exchange contribution, seen in Fig. 8(a). 

• The invisible decays of h, H, A and are expected, and plotted in Fig. 9(a). 
Future searches at the LHC, and in particular, searches at the ILC may reveal the 
true nature of the DM particle. 

For the co- annihilation scenarios as in II- A and II-B, 

• Although the "well-tempered" scenario with large Higgsino and Wino fractions is 
disfavored by the XENON- 100 data, the co-annihilation may still be a valid solution 
to obtain the correct relic density. There may be several light SUSY particles such as 
neutralinos, charginos, or stau, leading to many rich phenomena that can be explored 
at collider experiments. 

• For highly degenerate NLSP, NNLSP, the decays of H, A and as shown in 
Figs. 9(b)-(d) could lead to large invisible modes, making the collider search for DM 
very interesting. 

We conclude that understanding the nature of DM requires us to consider results 
from a number of different experiments. Future collider searches and the next generation 
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of direct detection experiments will likely cover the conventional parameter range of the 
MSSM if the LSP constitutes all of the DM. The recent exciting discovery of the SM-like 
Higgs boson, and searches for beyond the SM physics at the energy frontier will serve as a 
new "lamp post" and guide in DM searches complementary to what may be obtained from 
direct detection and astro-particle observations at the cosmic frontier. 
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Appendix 



A. Relic density calculation 

When the Hubble expansion H = a/a became much larger than the interaction rate V = 
n x (a" a v), the WIMPs (x°), once in thermal equilibrium with the rest of the Universe, 
decoupled from equilibrium. The number density of WIMPs at a time t, n x (t), is obtained 
by solving the Boltzmann equation 

1 d (a 3 n x ) 

-3 dt = ~M K - <] > (A-l) 



a 



where {a a v) is the WIMP annihilation rate averaged over velocities, and n eq is the equilib- 
rium number density of WIMPs: 



g measures the number of relativistic degrees of freedom, and T is the temperature. 

Define the dimensionless variables Y = n x /s and x = m x /T, where s is the entropy 
density given by 



45 

Here and henceforth, we adopt the natural units kb = h = c = 1. g s is different from g only 
at late times, after neutrinos have decoupled from equilibrium, and annihilation leads 
to the photons being heated relative to the neutrinos. However, dY/olx is very small at late 
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times, and with good accuracy, we may set g = g s when computing the relic abundance. 
The entropy per comoving volume is conserved, and hence 



d(a 3 s) _ d (ga 3 T 3 ) 
dt ~ dt 

We may then rewrite Eq. (A.l) in terms of Y as: 







dY 
~dt 



-(a a v)s[Y 2 -Y 2 



eqj i 



where Y e 



eq 



n 



eq 



/ s. From Eq. (A. 4), we see that ga 3 T 3 = constant, and therefore, 



-l = H + <L. 

T 3g 

The Hubble parameter H{T) is given by the expression 



(A.4) 



(A.5) 



(A.6) 



H 



8vrG 



I 1/2 



Pi 



8n 3 G^ 1/2 



90 



g l/2 T 2 



(A.7) 



Differentiating x = m x /T with respect to time, we find 



-T/T)x = Hx 
« Hx 



1 + 



1 



9s 
3Hg s 
1 d(lng 



3 d(lnT) 



(Ai 



where we simplified the second term on line 1 by substituting x ~ Hx (provided g s <C 
3Hg s ), and therefore g ~ —HT(dg/dT) [106]. Note that 5 changes significantly at the 
epoch of quarkdiadron transition. We may now rewrite Eq. (A.5) in terms of x: 



dY 

dx 



-(a a v)s(x) 



H{x)x 



1 



1 d(\ng) 
3 d(lnT) 



7T 

45G 



(a a v)m x 



9 



1/2 



1 1 eq 



1 



1 d(lng) 



3 d(lnT) 



(A.9) 



We solve Eq. (A.9) numerically to obtain the present day value loj once the form of 
g(T) and (a a v) are known. The dark matter relic density is then computed as: 



a v h 2 



m x Y s 
(Pcrit/h 2 ) ' 



(A.10) 



where so ~ 2893 cm 3 is the present day entropy density and p cr it ~ 1.054 x 10^ 5 /i 2 
GeV/cm 3 is the critical density. 

After performing a numerical integration of the Boltzmann equation as formulated in 
Eq. (A.9), we show the WIMP number density in Fig. 10(a) and the WIMP relic (mass) 
density in Fig. 10(b), for various WIMP mass values. The dark straight-falling line gives 
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Figure 10: Evolution with temperature and yielding the correct WIMP relic density Vl x h 2 = 0.11 
with illustrative values of the WIMP mass m x = 5, 100, and 1000 GeV, (a) WIMP number density, 
and (b) WIMP mass density. The equilibrium lines are for m x = 100 GeV. 

the densities if the particle keeps in thermal equilibrium with the environment for m x = 100 
GeV. It is known that the freeze-out temperature for a relatively light WIMP particle is 

xt = m x /T ps 20. (A.ll) 

The horizontal curves in Figs. 10(a) and (b) present the WIMP number density and mass 
density after freeze-out for m x = 5 — 1000 GeV, leading to the correct relic density. 
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